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‘Event Horizon, a null hypersurface defining the boundary of the black hole region of a 
spacetime, is not particularly useful for evolving black holes since it is non-local in time. 
Instead, one uses the more tangible concept of Apparent Horizon for dynamical black holes 
out there in the sky that do all sorts of things: evolve, merge and feed on the environment. 
Event Horizon, being a gauge-independent, global property of the total spacetime is easy to 
define and locate in the stationary case; on the other hand, Apparent Horizon depends on 
the embedding of the surface in spacetime and hence it is somewhat tricky to define. But 
for numerical simulations in General Relativity, locating the Apparent Horizon helps one to 
excise the black hole region and the singularity to have a stable computation. Moreover, for 
stationary solutions the two horizons match. Here we give a detailed pedagogical exposition 
of the subject and work out the non-trivial case of a slowly moving and spinning black hole. 


I. INTRODUCTION 


Stationary (Kerr) and static (Schwarzschild) black hole solutions of General Relativity 
have rather dull lives: stationary ones do the same thing, static ones do nothing as observed 
by an observer outside the black hole. While these vacuum solutions obtained in an isolated 
universe serve as our starting point for a more physical and detailed understanding of actual 
astrophysical black holes, the latter are almost never isolated: the black holes out in the sky 
have accretion disks, companion stars, neutron stars or black holes. Black holes feed on their 
environment and grow; in fact they are the most dynamical parts of the vacuum. As the 
first LIGO/VIRGO gravitational wave detection showed [I], black holes can grow feeding on 
other black holes: cannibalistic behavior of these objects-highly curved vacua-could explain 
the existence of intermediate mass black holes. 

As the astrophysical black holes evolve, concepts such as the Event Horizon defined easily 
for eternal black holes are not clearly adequate for us, the transient observers. Recall that 
the Event Horizon (H) of a stationary black hole is a co-dimension one null hypersurface 
in the totality of the spacetime defined as the boundary of the black hole region from which 
time-like or light-like geodesics cannot reach future null infinity P]. Stated in a different way: 
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it is the boundary of the region which is not in the causal past of the future null-infinity. 
This says that the Event Horizon is a global property of the totality of events which is all 
of the spacetime. Therefore, one cannot locate the Event Horizon with local experiments in 
a finite interval of time. In this respect, it is apt to say that the Event Horizon Telescope 
detected the environment of the black hole from which one can see at best the cross-section 
of the Event Horizon, not the Event Horizon itself. 

For dynamical black holes one invents the more useful concept of the ‘‘Apparent Horizon” 
[3], a co-dimension two spatial surface (hence local in time), which, unfortunately, in general 
does not carry geometric invariant data as the Event Horizon but it contains sufficient 
information regarding the possible formation of an Event Horizon in the future that it pays 
to describe it in detail. In numerical relativity computations, detection of a black hole region 
is best done with Apparent Horizons. Within the context of General Relativity, existence 
of an Apparent Horizon implies the appearance of a future Event Horizon outside of it. 
Therefore, one can excise the region inside the Apparent Horizon (that also includes the 
singularity) for the stability of the computation since nothing will come out of that region 
in classical physics. For modified gravity theories, an Apparent Horizon need not be inside 
the Event Horizon (See the discussion and references in [4]). 

Our task in this work is to give a detailed definition of the Apparent Horizon and some 
related concepts and apply it to slowly rotating and moving initial data which was recently 
given in [5]. The layout of this work is follow: in section II we introduce the necessary 
tools for the defining equation of an Apparent Horizon as a co-dimension two spatial hy- 
persurface in n dimensions and use the ADM decomposition of the metric to arrive at an 
equation in local coordinates, in section III we consider a conformally flat initial data for 
n = 1 + 3 dimensions for which the momentum constraints can be solved exactly following 
the Bowen-York construction [6]; and we solve the Hamiltonian constraint for slowly mov- 
ing and spinning initial data and compute the properties of the Apparent Horizon. In the 
Appendix we expound upon some technical points alluded to in the text. 


II. DERIVATION OF THE APPARENT HORIZON EQUATION 


As stated above, the Event Horizon of a black hole, as a null hypersurface, cannot be 
determined locally: one has to know the total spacetime to define it. On the other hand, the 
Apparent Horizon can be determined locally in time. For this purpose, we need to define a 
congruence of null geodesics and its expansion. Our notations will be similar to those of the 
excellent lecture notes [2] [7]. 

As shown in Figure 1, we have an n dimensional spacetime manifold .@, with a co- 
dimension one spacelike hypersurface X, that is dim = n — 1; and we introduce a co- 
dimension two subspace S, dim S = n — 2. Let n” be a timelike unit vector orthogonal to 
yi: 

a iy = —1, (1) 


and s” be a spacelike unit vector orthogonal to S 
ss, =1. (2) 
We impose the condition that n and s-vectors are perpendicular to each other 


ns, = 0. (3) 


Figure 1: On an n dimensional spacetime manifold .@, we introduce a co-dimension one spacelike 
hypersurface © and also a co-dimension two subspace S. Timelike unit vector n” is orthogonal to 
X and the spacelike unit vector s” is orthogonal to S. k” and Æ denote the ingoing and outgoing 
null vectors respectively. 


Instead of these two vectors, one can also work with the ingoing null vector k” and the 
outgoing null vector 4, defined respectively as follows (see Figure 2) 


k” := 3 (n# — s”), L = nt H F. (4) 


The induced metric on the hypersurface X is 
Yuv = Juv F Nyy, (5) 
while the induced metric on the subspace S reads 
Quv = Yuv — SuSv = Juv + Nw — SySv, (6) 
where u,v run over the spacetime directions. The important concept here is the extrinsic 
curvature of both of these surfaces. For the hypersurface X, we have 


Kw = —YneueV Nn’, (7) 


where V,, denotes the covariant derivative compatible with the spacetime metric, Vg.) = 0. 
From a more geometric vantage point, our definition is as follows: given two vectors (X,Y) 
on the tangent space at the point p, that is T,X, and n being the unit normal to X, then 
the extrinsic curvature of X is defined as K(X,Y) := —y(Vxn,Y). So in local coordinates, 
one can take X = 0,, Y = 0, to get Kuw := K(0,,0,) = —7(Va,n, 0.) which matches (7). 
The minus sign is a convention. Equivalently, one has? 


Kye = V pno = yl Vay: (8) 
Similarly, we define the extrinsic curvature of the (n — 2)-dimensional space S as 


kuv = —quolvp VT P, (9) 


2 One can also define the acceleration as a := Vpn. 


Figure 2: The unit vectors n“, s“ together with the ingoing null vector k“ and outgoing null vector 
é¥ are shown. 


and using the definition of the induced metric (6) one obtains 
Ke = = V phe — NIV pSo = TA Viste — Tata VS + SpSV oSv + SpNuSoNp & (10) 
One defines the expansion of the out-going null geodesic congruence as 
Oig = dV ilu, (11) 


which is the divergence of the null geodesic congruence along its propagation in the outgoing 
null direction. Using and the extrinsic curvatures of the hypersurface and the surface, 
we can recast the expansion of the null geodesic congruence as 


Oa =K+k+ (nn — s"s”) (Ky + kuw). (12) 


Since nn” Ky, = 0 = n”n” kyu = s” kuv, Oq reduces to the following neat equation as 


Ov) =—-K—-—k+ Bee K pv- (13) 
Equivalently one has 
Ove) = =q” (Kw + k); (14) 
or 
Ow =q” (Ki + ka), (15) 


where the 7,7 indices run over coordinates on the hypersurface X. 

The expansion ©(y) is employed to define the very important concept of a trapped surface. 
An outer trapped surface on © is a closed (that is compact without a boundary) co- 
dimension two surface such that for outgoing null geodesics orthogonal to the surface, one 
has Ow) < 0 everywhere on the surface. The subset of X that contains the trapped surfaces 
is called the trapped region T, a co-dimension one surface. Finally, Apparent Horizon is 
the boundary of the trapped region (an obviously spatial surface) which we shall denote by 
S := OT. By definition Apparent Horizon is a marginally outer trapped surface (MOTS) 
and satisfies the Apparent Horizon equation: 


Ove) =—-K—k+ ae K w = 0. (16) 


It is clear that for the case of time-symmetric initial data (A, = 0), the Apparent Horizon 
becomes a minimal surface since k = 0. 

Now that we have defined the Apparent Horizon, given a metric in some coordinates, to 
proceed we need to lay out in detail how is expressed in terms of the metric functions. 
For this purpose we choose the ADM decomposition of the metric [8]. 

Let N = N (t, 2") be the lapse function and N’ = N‘(t,x’) be the shift vector, then the 
line-element reads 


ds? = (N,N* — N*)dt? + 2N,dt dx’ + y;dz' dz’, (17) 

or in components one has 
go = N;N* — N?, Joi = Ni, Jij = Tijs (18) 

with the inverses given as 
g? = —N-~, g” = NİN, gi = 4” — N'NIN-?. (19) 

Using the definition (5), one has 

Jij = Vig — MUNG = Vaz, (20) 
hence n; = 0. Similarly the relation g = 77 — n'nd = g — NNI N? yields nê = +N'/N. 
Since n” is a timelike vector, using n; = 0, one has no = +N and choosing for N > 0, we 


choose the plus sign for the future-directed time-like vector to arrive at 


në = & -7) ; Ny = (-N, 0). 


We can work out the additional relations between the spacetime metric g and the metric of 
the hypersurface y as follows 


(21) 


goo = Yoo — Nono = NN‘ — N°, (22) 
which yields yo9 = N;N*. And similarly 
Joi = Yor — Non; = N; (23) 


yields yo; = N;; and from the inverse metric relations, one obtains 7°" = 0. 

Similar computations for the co-dimension 2 spatial subspace S, after using the condition 
ns" = 0, yield s° = 0 and q™ = 0. 

Now we go back to and express it for the apparent horizon as 


From (8), one has Kj; = —Vinj = —NTY, and we obtain’ 


q (Kij = OiSj + Ty.s0 + Ti sn) = 0. (25) 


3 Note that, explicitly, we have Ki; = sw (DiN; + D;N; = tViz). 


We denote the Christoffel connection of the induced metric y as *T%,. Then substituting the 
corresponding components of the Christoffel connection one has 


k _ Spk Ne 
and we arrive at 1 
qg” (Ka = Dis; + wisi (N's. a +0) = 0, (27) 


where D; denotes the covariant derivative compatible with the spatial metric y, DiYjk = 0. 
One has N*s; — so = 0, so then the equation defining the apparent horizon becomes 


Before we start working out an example, let us note that there is another simple expression 
of the 0g in (11) and hence equation (28). One can show that (see section A of the Appendix 
for the proof) Lie-dragging the metric on S along £ yields exactly the expansion: namely, 
one has 


1 
Ae) = GV uly = gf Lelu, (29) 


where £e denotes the Lie-derivative along the vector £. In section B of the Appendix, f(g is 
derived from the minimization of the area along the outing null direction which also yields 
a complementary physical picture. 


II. APPARENT HORIZON DETECTION 
A. The equation in explicit form 


From now on we shall work in n = 1 + 3 dimensions. Assume now that the local coordi- 
nates on È are denoted as (r, 0, ) and that the location of the Apparent Horizon depends 
both on 0 and ¢. The equation to be solved is 


q” (As; — Dos, — Kij) =). (30) 
Assume that the surface S can be parameterized as a level set such that 

O(r,0,¢) :=r—h(6,¢) =0, (31) 
with h being a sufficiently differentiable function of its arguments. Since sê is normal to 


the surface, one has s; ~ 0;®; and because it is a normal vector, let s; := Að; and define 
mi := 0;®, which yields 


me a(i, E -ðh ). (32) 
To proceed further, let us take the metric on © to be conformally flat as in [6] 
1 0 0 
gaa 07 O |, (33) 


0 0 r*sin?@ 


then one has 
s' = Ce =% Əh, —¥** Ash) , (34) 


with the normalization factor given as 


-1/2 
X= (a +A? (ðh)? + 7°*(Bph)?) , (35) 


As should be clear at this stage, the Apparent Horizon equation will be a rather complicated 
non-linear partial differential equation with little hope to yield an exact analytical solution. 
Let us further assume (following [6]) 77 K; = K = 0, which is called the maximal slicing 
gauge. Then (30) reads more explicitly as 

yom; — Teme — \*m'm!0;m; + Ammim, + Am'm! K; = 0. (36) 


After working out each piece in a somewhat tedious computation, one arrives at 
1 
= Oph — 4 gh — 5 (08 — ey a = at doh? dae) 


+42 (0%) aoh) ah + (79°)2(Osh)2O2h + 2yH%-y"dphOohdodsh 
y2 
F (Ga + (4) (Ogh)?OnV00 F (P (Oh) Yo 


2 
- (8h) Bohl)? Aerie) 
HAÇ Kr + (POA) Koo + (P (Ooh)? Kop — 27" Doh Ke 
Ly" Osh K re + 290 y*Aghdoh Kos) = 0. (37) 


Given the metric y;; and the extrinsic curvature Kj, one can find numerical solutions of 
this equation up to the desired accuracy. Our goal here is to find approximate analytical 
solutions to some physically reasonable initial data which must satisfy the Hamiltonian and 
the momentum constraints on the hypersurface X which we discuss next. 


B. 1+3 form of Einstein equations 


The Hamiltonian and the momentum constraints on the hypersurface X follow from 
Einstein’s equations as 


—"R— K? + Ky K” — 2«Ton = 0, 
2D, KE — 2D;K — 2kTri = 0. (38) 
We chosen K = 0 and consider the vacuum case with Tu» = 0. Of course this initial data 


evolves in time and the remaining parts of the Einstein equations written as a dynamical 
system are given as 


LoT 
at —_ 
ð 
Ka = N (Rij — "Rij — KKy + 2K KË?) + ZgKiy — DDN, (40) 


where -Y ș is the Lie derivative along the shift vector N’. Derivation of these well-known 
equations can be found in many textbooks, see our derivation in [9]. 


IN Ki; + DiN; + DjNi, (39) 


C. Conformally flat Bowen-York type data 


For a conformally flat hypersurface © (yj = Yt fi; with f being the flat metric in some 
coordinates), the constraint equations (together with the "maximal slicing" condition 
K = 0) reduce to a non-linear elliptic equation and an easily solvable linear equation, 
respectively given as 


he, thes 1 X 
D;D'y = -3V Kip (41) 
D'R; = 0, (42) 


Bowen and York [6] gave the following 7-parameter (p;, a, Ji) solution to (42) on R? whose 
origin is removed: 


3 3a? 
i= 52 (vin; + pjni + (ning — fiz)p n) Fia (vin; + pjni + (fig — Sninj)p - n) 


3 
a ki (erun; + enti) (43) 


A 


where r > 0 is the radial coordinate, n’ is the unit normal on a sphere of radius r (not 
related to the unit normal to ©); € = +1 and p- n = p*n x. At this stage, one should note 
that the physical meaning of the parameters (p;, a, J;) is not clear; secondly, linearity of 
means that each bracketed term solves the equation separately. For the sake of simplicity, 
we shall choose a = 0 in what follows. 
Here we follow [5]. We shall need the following expression for the right-hand side of 
os one 9 18; = 18 ;> > 
Pena (P +28- 0)) + = (J xñ) B+ T (J xñ) (Tx a). (44) 
Inserting this expression to (41), one arrives at the complicated Hamiltonian constraint 
which can only be solved exactly after making several assumptions. We shall not go into 
that discussion which was given in [10] in some detail. 


D. Conserved quantities 


To understand the physical meaning of the parameters in the solution, we shall assume 
that the spacetime is asymptotically flat, hence the conformal factor behaves as 


w(r) =14+O0(1/r), as r —> OO. (45) 


Then one has the conserved total momentum associated to © easily written as a boundary 
integral on a sphere at spatial infinity: 


1 . 1 s 
8T S2 8T S2 
Observe, from the second equality, that only the leading term in the conformal factor is rel- 
evant for this and the following computation. The total conserved total angular momentum 
can also be found easily as 
1 


1 7 - A 
= grn tik I. dS nisi K” = gn i. dS n ri K". (47) 


Given and (47), it is straightforward to compute the integrals for the extrinsic curvature 
(43) which at the end yield P; = p; and J; = Ji. So for the computation of these two 
quantities, let us note once again that, the full form of the conformal factor is not needed; 
one only needs to know its behavior at infinity, that is the O(1) term. 

From these two conserved quantities, one can see that physically the assumed extrinsic 
curvature belongs to a a self-gravitating system (a curved vacuum) with non-zero mo- 
mentum and angular momentum. To compute the total mass-energy, the ADM energy, of 
the system, the O(1) term of the conformal factor is not sufficient. For that computation 
we keep the next order term and assume 


w(r) =1+ 5. +O(1/r?) as r 00. (48) 


Then defining the deviation from the background as hij := (Y4 — 1)di;, the ADM energy 
simplifies as 


Esp = 7 f asin; (Bh? — dhi) = -z> | dS ni du, (49) 
SZ S2 


whose explicit evaluation for yields E4pm = E, which of course at this stage is almost 
a tautology: we have to find the constant Æ by solving the Hamiltonian constraint. 


E. Approximate solution of the Hamiltonian constraint for a boosted slowly 
rotating gravitating system 


To solve the elliptic equation (41) using (44), let us take k to be the direction of the 
conserved angular momentum and choose p to be in the xz plane (this is just a choice of 
the orientation of the coordinates and no generality is lost) 


> 


J = Jk, B= psin O91 + p cos Ok, (50) 


with 6o a fixed, conserved angle. Then the Hamiltonian constraint becomes 


A Ai 9J QJ? 9p2 
DD =4 (2a sin 6 sin ¢@ — qa sin’ 6 — ai 


(1 + 2(c, sin 6 cos @ + c2 cos a) ; 
(51) 
where cı := sin fo, C2 := cos 0o. 

Needless to say, an exact solution of this equation is hopeless, therefore we shall search 
for the lowest order perturbative solution assuming an expansion in terms of the momentum 
and spin which corresponds to a curved 3-surface with a small linear and small angular 
momentum. In [IJ] the slowly spinning case with no linear momentum was solved in the 
leading order; and in [12] slowly moving without spin was solved and in [5], both motions 
were considered at the leading order. We now present this solution. 

A cursory inspection on the right-hand side suggests that one should have a double series 
of the form 


plr, 0, p) = PO + Py + pP + Jp + Ol, Jt PI), (52) 


10 


where the functions on the right-hand side depend on (r,6,@). At the zeroth order, one has 
the usual Laplace equation oe 
D;D'y© = 0, (53) 


which needs boundary conditions to be uniquely solved. The following boundary conditions 
as employed by [I2] are apt for the problem at hand: at spatial infinity one demands 


lim v(r)=1, yr) > 0 (54) 
and near the origin one has 
lim, y(r) = ¥, (55) 


where ~© might have a singularity at the origin. In fact the zeroth order solution satisfying 
these boundary conditions reads 


yO =1+5, (56) 
r 


where a is a constant of integration at this stage not to be confused with the one in (43). 
The equations at the next order are 


D,Dip = r: sin? 0 pay (57) 

D; Dip?) = a sin ô sin pooo (58) 
‘ 4 (r+a)” 

D, Dip = a (1 + 2(cı sin 0 cos @ + cg cos ay”) ee (59) 
‘ 16 (r+a)’ 


These are linear equations whose solutions can be found with the help of the following 
spherical harmonics : 


1 3 | 5 
Yo (0, ¢) = Vir YP (0, ¢) E JŽ os 0, Yp (0,9) E Ten cos’ 6 a 1), 


1 3 
Yr (0, $) = \/ sin Asin 6 Y3(0,¢) = I sing cos cosg, Y (0, o) = y g Sin A cosg. 


Then a close inspection of suggests that the proper ansatz for w/) should be of the 
form 
VO (r,0,6) = YD O, p) + VYRO, 4), 
from which the solution obeying the boundary conditions can be found to be 
(a* + 5a°r + 10a?r? + 5ar? + r*) r? 
40a3(a + r)> A0a(a + r)5 


wr, 6,6) = (3 cos? 6 — 1). (60) 


To solve one should take 
P(r, 0, p) = Vy (r)¥o (0, p) + VEP (YTO, 4), 
for which the solution obeying the boundary conditions is 


cır (a? + 5ar + 10r?) 
80a(a + r)> 


YUP (r 0,6) = sin @ sin @. (61) 
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The Y”) equation is similar albeit slightly more complicated: the proper ansatz reads 


Y = p (r)YD (0, P) + YP YEO, 0)? + DP (YO, 0) + BS YRO, 6), 
from which four equations follow whose solutions are as follows: 


Jr (84a? + 378a°r + 653afr? + 514a3r? + 142a?r4 — 35ar° — 25r®) 


(p) es 
o (r) 80ar?(a + r)’ 
21ra a 
= ] 62 
20r’ °8 a+r’ (62) 
and 
Te (84a° + 378afr + 658a°r? + 539a?r3 + 192art + 15r® 
pl?) (r) _ 1 
g 40r?(a + r)* 
21nac? a 

: 63 
lore $a ee 


yP (r) can be obtained from with the replacement cj + \/ cic. and we ) (r) can be 
obtained from (63) with the replacement c? > c3 . All these pieces can be combined to 
get w) at this stage, but a depiction of the final result is redundant since all the parts are 
given above and the final expression is cumbersome. We have now all the information at 
our disposal to compute the relevant quantities defined on © including the location of the 
Apparent Horizon. 

First let us revisit the ADM energy computation which we started above: We need the 
dominant terms up to and including O(2) in W(r,0,¢). A quick power series expansion 
yields 


a J? 5p? 1 
=1 
ule) ss r a 40a3r  32ar OC): (pa) 


in which the Jp term appears at O(4) and therefore makes no contribution to the energy. 
Then from (48), the ADM energy of the solution follows as 


r2 


J? 5p 
ee. (65) 


Pe 
ADM = AAt a pe 


So one can immediately see that for vanishing spin and vanishing linear momentum (that is 
the case of the Schwarschild black hole written in the isotropic coordinates) the constant a 
is related to the mass of the Schwarschild black hole mass as a = M/2. 


F. Apparent Horizion area and the irreducible mass 


While studying the efficient processes of extracting energy from rotating black holes, 
Christodoulou [13] realized* that there is an irreducible mass Mirr which is related to the 


4 Note that after Hawking’s area theorem [14] which came later than Christodoulou’s observation, it became 
clear that there must be an irreducible mass at the classical level. 
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area Apy of a section of the event horizon via 


| AnH 
Mir := 4] ——. 66 
167 ( ) 


For a moving, rotating black hole, the total energy was obtained in [13] as 
2 


4M?..’ 


irr 


EP = Mgr +p + 


(67) 


in which the physical meaning of each part is clear. 

Since we have a dynamical, evolving system, we have at our disposal the area of the 
Apparent Horizon, not a section of the Event Horizon. But, following [I2], a good approxi- 
mation to Mir can be given with the help of the area of the Apparent Horizon via 


_ {Aas 
Mie = 4/ Te (68) 


As we shall see, this definition yields the correct expression for the energy of our system 
obtained from an expansion of (67). But first we need to find the location of the Apparent 
Horizon, namely solve up to the accuracy we have been working with. That area is 
given simply as 


2T T 


À= / dg | A E (69) 
0 


0 
which yields the following exact form: 


20 T 


, 1 2 1 A 
Aan = fofa sin 0 yt h? (1 + $5 (dh) + Fesin2 0 (Osh) ) (70) 
0 


0 


Hence to get the area, all we need is to find the location of the Apparent Horizon up to first 
order in the spin and momentum. This suggests the following ansatz: 


h(6,¢) = h? + ph? + Sh? + O(p?, J’, Jp), (71) 


where 
ðh = 0, ðh? = 0 = h? = 05h". (72) 


Ignoring the terms such as (Ogh)?, (Osh)? and ðshðsh, (37) reduces to 


1 sor, bs 3p 
Ogh + — 55 Osh + cot 03h — 2r — dr m TE oh- ya 


At the zeroth order, O(p?, J°), it yields 


(c sin 0 cos @+ C2 cos 6) =). (73) 


Ont 


1+ 2r 
Y 


0, (74) 
with Y = 1 + £ ; setting r = h, one finds 


h? =a. (75) 
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This solution identifies the parameter a as the location of the apparent horizon at the lowest, 
dominant, order. For example, for the Schwarzschild black hole h = 2M (as noted above) 
would be the exact solution for which the Event Horizon and the Apparent horizon coincide 
in these conformally flat, isotropic coordinates. 

At O(p) we have an inhomogeneous, linear Helmholtz equation on a sphere (S°), 


GE ey + cot 609 — 1) h? = Ž (a sin cos @ + c2 cos 0), (76) 


while at O(.J), we have a homogeneous one: 


(+ = srg + cot bay — 1) W = 0. (77) 
Therefore, we have to find everywhere finite solutions of the following equation 
(Vi. + k) £ (0,6) = g (0,9), (78) 
where Va is the Laplacian on S$? given as 
V2 i= Og + cot 809 + — z ag% (79) 


One can employ the Green’s function technique to solve this problem. For the Helmholtz 
operator on the sphere, the Green function G(#, #’) is defined as 


(Viz + MAF 1)) G8, 2’) = 6 (@ — 4), (80) 
which can be found as an infinite series expansion (for example, see [15]) 


Bs i i & 1 Peers aah ieee 
Gaz) = ing 2 aly T'(-A) ie +1) ( 9 ; 


(81) 


where ĉ = sin 0 cos di + sin 0 sin ġj + cos 0k and ĉ’ is a similar expression with some other 


0 and @. Employing this Green’s function with À = as one finds the first non-trivial 
correction to the location of the Apparent Horizon as 
1 
h? = —— (c sinf cos ọ + cg cos 0), (82) 
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and h? = 0. Therefore the apparent horizon is perturbed from the zeroth order expansion 
to 


r= h(0,ġ)=a— E (singo sin 0 cos @ + cos ĝo cos 0), (83) 


where, recall that, 0o is the angle between the linear momentum and the spin vectors. So the 
magnitude of the spin vector does is irrelevant at this order for the location of the Apparent 
Horizon, but the angle it makes with the momentum vector is relevant. In Figure 3, we 
plotted and example of how the shape of the horizon looks like. There is a dimple on the 
sphere whose size depends on the ratio p/a which we took to be large to see the dimple. In 
the limit 69 = 0, h reduces to the form given in [12]. 
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Figure 3: Shape of the Apparent Horizon when the angle between p and J is 45 degrees; and 
p/a = 8/2 which is outside the validity of the approximation we have worked with. 


Let us now evaluate the area of the Apparent Horizon from which at the end yields 


llr J? 
Aan = 64ra? + 4p? + Ba? (84) 
Thus the irreducible mass Mirr turns out to be 
2 2 
p 11J 
Mirr =2 Gum DaN 3) 85 
a+ Iba © 320a? (85) 
so comparing with the energy, E4pj, we have 
2 2 
p J 
Bapm = Mir + IM. + VE (86) 


which matches the result of Christodoulou at this order. 


IV. CONCLUSIONS 


We have presented a step-by-step construction of the Apparent Horizon equation which 
is of extreme importance in black hole physics; and described in detail how it correctly 
yields the expected results, such as the irreducible mass, for a slowly moving and spinning 
black hole. For stationary black holes the event Horizon and the Apparent Horizon coincide. 
This exposition is of a pedagogical nature with details given in the Appendix including the 
derivation of the null Raychaudhuri equation which we have not used in the text, but added 
for more insight for the expansion of a null geodesic. We have skipped some interesting 
issues such as: numerically solving the case with no symmetry; multi black hole initial data; 
the proof that when the dominant energy condition is satisfied, the topology of the Apparent 
Horizon is that of S?. For other nice expositions regarding horizons and related concepts 


see [7] 16] 7]. 
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V. APPENDIX 
A. An equivalent definition of the expansion Oç) 


Here we give a proof of the second equality in (29): we have 
1 
Ov) = GV ky = g1 Lelu, (87) 


where the first equality is identical to the definition of the O(g. Starting from q” Lequv, one 
can easily arrive at the expansion Og). The construction is as follows: 


q” Lequv = CV oquv Ve + qouV l). (88) 


The first term on the right hand side automatically vanishes. To be able to see this explicitly, 
we express the metric quv in terms of the spacetime metric guv 


OY oquv = 40V o (guv + kulo ey) « (89) 
Since Voguv = 0, the non vanishing terms are 
EN oquv = 4C (Ku Voly + 6V oku + ky Volu + fuVoky) ; (90) 
where q” k, = 0 = q” l„, and so one gets 
q UNV oduy = 0. (91) 


Now let us evaluate the second and third terms in (88) (which contribute equally). We can 
write 
q” Ge Vue = q” (Gan + koly + kulo) Vl. (92) 


Using q” k„ = 0 = q™” l, again, the last expression reduces to the following 


Q” qovV pkl = Jor) ul =a VN pho. (93) 
Then becomes 
q” Loduy = 24” V plv, (94) 
and one ends up with 
1 
g1 Ledu = aN a = Oe), (95) 


which is the expression we wanted to prove. 
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B. Derivative of the area along the null vector field 4# 


To gain a better physical insight to the expansion of the null geodesic congruence, let us 
show that when one takes the derivative of the area of the cross section along the null vector 
field 4”, the expansion ©) can be directly obtained as the integrand. On the surface S, let 
us start with the area formula 


A= J dS /q, (96) 


of which the derivative along £ yields 


1 
HO,A = l dS ai= 5 i dS Gli qn dav: (97) 


Equivalently one can express the result in terms of the Lie derivative along the vector field 
é using 
g Lidar = g l Oaa + 29 qual”. (98) 


Since the null vectors @“ and k” are the elements of the compliment of the subspace S, one 
has k* = 0 = @*. By definition (6) we obtain 


Qn = On + Rye + RE, = On, (99) 


and similarly 
qr = g ere eee = g”. (100) 


So that ae ey = Ôn; and the last term in becomes 
gq” qurOal" = 820, = Age? = 0. (101) 


Then gL qap reduces to 
q” Leda = ql Ondab; (102) 


which can be related to "0,,/q via 
Poya = 5 Ja Lega (103) 
Now we can rewrite in terms of Lie derivative 
£0,A = fas Vd ; a” Ldap: (104) 
In order to show the appearance of the expansion ©) explicitly, we should use the spacetime 


coordinates, recall that we have Oq) = q“”Leqy/2, instead of the coordinates on the co- 
dimension two surface S. It is straightforward to write 


ga Ledar = SOL edar = Lelu- (105) 


Multiplying this with qfq% one obtains 


aed Leda = a Leda, (106) 


which in terms of Kronecker delta functions reads 
gg 528, Legas = aq Lequy, 
and yields 
q9 Legan = WG Leu: 
We multiply the last identity with q°°. Then we find the identity 
geio 7 Leda = 9 akg Leduy 


and so one arrives at 
ab — {0p jb AV 
q” Legas = 1° OG dp Led: 


where q°?qhq, = q"”. Finally we end up with 


q” Legas = 9” Ladue 
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(107) 


(108) 


(109) 


(110) 


(111) 


This proves that the expansion Op appears in the change of the area along the vector field 


4”. The final expression is therefore 


1 1 
L0, A = fas Vd 5 q” Lequy = fas vaq 5 oe Lits= fas va O- 


So setting Op = 0 to define the Apparent Horizon boils down to setting (“0,,A = 0. 


C. Null Raychaudhuri equation 


(112) 


The form of already suggests that one can define a tensor whose trace is the ex- 
pansion. Here we explore this tensor and obtain an expression for the change of the null 
expansion along the null direction £ as well as the null Raychaudhuri equation [2]. So let us 


introduce the deformation tensor O, as 


1 o 
Ow = Zh eop; 


such that 
Ow) = g” Ow = q” NV lo. 


Carrying out the Lie-derivative in (113), one has 
Opm = V ple = yl a Volby, 


with 
Wp == Vile — Kk Whale 


which is called the rotation one form. 
In what follows, we will make use of the Ricci identity 


VaV — VV = Bil. 


From (114), one has 
Oe) = Vie" + KLEN ee 


(113) 


(114) 


(115) 


(116) 


aa 


(118) 
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Here we assume that £ is a geodesic null vector but not necessarily affinely parameterized 
so that 


LEV hy = Key, (119) 
where « is a function on spacetime. Using k’¢, = —1, one has 
Ow = Vie" — K. (120) 


So we have the following two equations: 
VÆ = Ow +k, Valu = Op Foil — lak V obv. (121) 
Substituting these in (117), one has 
Vu (O, t + wlt — keV ob") — Val O +6) = Ril, (122) 
which more explicitly becomes 
V pO, + LAV pwy HUV lt — keV LV ly — OV (ke Vol") — Vi(Ow +6) = Bil. (123) 


We use the expressions (121) one more time and reexpress the third and the fourth terms 
to obtain 


VO,” + EV pwy + w,(Qw + &) — kV (Quy + wylo — luk’ V yr) 
L, V (k7 Vol") — V (Op +k) = Rab. (124) 


Since l Vol” = 0, the last term on the first line automatically vanishes. Contracting the 
final expression with @” and using the fact that it is a null vector, one arrives at 


OV Oy" + LEY wy + u(y +k) — LO, kV ol! — LV (Ow + K) = Ryle’. (125) 


It is easy to show that the contraction of the null vector @” and the deformation tensor 
identically vanishes, ©,” = 0. Then one has 


PVO," = -0p 9, (126) 


and also 
UN pwy = 00 pK K’. (127) 


Inserting these expressions in (125) we get 
— Op O% + KO() aa LV OQ = Ral’. (128) 


The first term 0,0” can be written in terms of the trace free shear tensor ouy 


Ow = Ow = Av P(e) (129) 


n—2 


as follows 1 
O, O” = Tua” + Bo, (130) 
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where o,,0 = ono. Then (128) becomes 


1 
LY Oio = VG) = KO) = get = Taa” = m , 


The null vector field 4” is oriented in the future direction. Therefore the last equation, 
known as the null Raychaudhuri equation, is an evolution equation for the expansion Oy). 
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